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The synthetic organic compound A-(BETS)2FeCl4 undergoes successive transitions from an antifer- 
romagnetic insulator to a metal and then to a superconductor as a magnetic field is increased. We 
use a Hubbard-Kondo model to clarify the role of the Fe 3+ magnetic ions in these phase transitions. 
In the high-field regime, the magnetic field acting on the electron spins is compensated by the ex- 
change field H e due to the magnetic ions. This suggests that the field-induced superconducting state 
is the same as the zero-field superconducting state which occurs under pressure or when the Fe 3+ 
ions are replaced by non-magnetic Ga 3+ ions. We show how H e can be extracted from the observed 
splitting of the Shubnikov-de Haas frequencies. Furthermore, we use this method of extracting H e 
to predict the field range for field-induced superconductivity in other materials. 
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The discovery of magnetic-field-induced supercon- 
ductivity |l| in the two-dimensional compound A- 
(BETS) 2 FeCl 4 [where BETS is bis(ethylenedithio)- 
tetraselenafulvalene] is a new example of the rich phase 
diagrams of organic molecular crystals [||. Whereas 
previously pressure or chemical substitution has been 
used to change the electronic properties of these organic 
materials, it is remarkable that this compound under- 
goes successive electronic phase transitions as the mag- 
netic field is increased. Below a temperature of 8 K, 
A-(BETS)2FeCi4 is an antiferromagnetic (AF) insulator 
||. As a magnetic field is applied, it undergoes a first- 
order transition to a metal at II T. Close to this field, 
the magnetic moments associated with the spin 5/2 of the 
Fe 3+ ions undergo a transition to a polarized paramag- 
nct. If the magnetic field is parallel to the layers, there is 
a transition to a superconductor at 20 T H , which is then 
destroyed above 42 T pi . The magnetic ions are essential 
to this behavior, since the compound with non-magnetic 
ions, A-(BETS)2GaCi4, is, in contrast, a superconductor 
at zero field H, despite very similar crystal structures 

In this Letter we focus on three questions: (i) Why 
does the inclusion of magnetic ions change the ground 
state from a superconductor to an insulator? (ii) Is the 
magnetic-field-induced superconductivity due to the the 
Jaccarino-Peter effect jjj (as has recently been proposed 
M), where the external field is compensated by an in- 
ternal exchange field due to the magnetic ions? and (iii) 
Does the Jaccarino-Peter picture survive if one takes into 
account the spin fluctuations associated with the mag- 
netic ions? 

Recently, Ziman introduced a two-dimensional Hub- 
bard-Kondo model in order to understand question (i) 
p|. The model takes into account the four conduction 
bands associated with layers of BETS molecules (four 
HOMO orbitals per unit cell), a Kondo coupling be- 
tween the localized S = 5/2 spins and the conduction 



electrons, and the Coulomb repulsion between two elec- 
trons on the same BETS molecule. Ziman found that for 
small electron-electron repulsion the periodic potential 
due to the magnetic ordering (found self-consistently) at 
low temperature opens energy gaps on the Fermi surface 
||. A magnetic field, by aligning the moments, destroys 
the periodic potential, restoring the Fermi-surface. How- 
ever, to suppress the entire Fermi-surface, this needs a 
Kondo coupling, J > 6 meV, which is larger than the 
estimates that we extract from experiment below. More- 
over, the system seems to have quite a large electron- 
electron repulsion, as suggested by comparison with the 
k-(BEDT-TTF) 2 X family g. In this case, we show first 
that the system without the magnetic ions may be close 
to a Mott transition. Then, the Kondo coupling with 
the magnetic ions can drive the system into the insulat- 
ing phase in order to gain some magnetic energy. These 
two scenarios of the metal-insulator transition lead to 
different physical pictures (band-insulator versus Mott- 
insulator) . 

Question (ii) has to be carefully examined. Although, 
it is clear that the magnetic ions can in principle pro- 
duce an exchange field H e that can compensate the ex- 
ternal field, it is desirable to know the precise magnitude 
of H e . We show how to extract it from the observed 
magnetic oscillations M . This allows us to rule out alter- 
native proposals such as spin-triplet superconductivity, 
field-induced dimensional crossovers, or superconductiv- 
ity mediated by spin fluctuations in the local moments. 

Previous estimates of J involve considerable uncer- 
tainty. In the high-temperature metallic phase, the ex- 
change leads to an RKKY interaction between the lo- 
calized spins, Jo = J 2 x(Qaf) where xaf is the elec- 
tronic spin susceptibility at the wave-vector of the AF 
correlations. The high-temperature magnetic suscepti- 
bility gives an estimate of Jo ~ 0.2 meV pf. To obtain 
the coupling J from this approach, we need to know the 
electronic spin susceptibility x{Qaf)- Using the free- 



electron band structure, x(Qaf) = 80 (eV) _1 || gives 
\J\ = 1.5meV. Hotta and Fukuyama [Q suggested that 
the Kondo coupling comes from superexchange processes 
leading to an antiferromagnetic coupling (J > 0). They 
estimated J ~ 1 meV, using hopping integrals found from 
Hiickel calculations and assuming a value of 2 eV for the 
splitting between the d orbital of the Fe 3+ and HOMO 
orbitals. 

Mott insulator. We first argue that the materials with- 
out the magnetic ions are close to a metal-insulator tran- 
sition. From the experimental point of view, the effect of 
the anion in A-(BETS)2GaBr 2 Cl4_ z is to drive the elec- 
tronic system from a superconductor for z < 0.8 to an in- 
sulator z > 0.8 pifl . As the crystal structure is very sim- 
ilar in both cases, this means that a small change in the 
electronic parameters (estimated to be smaller than 5% 
JL2JI) yields two different phases. Hence, the electronic 
system without magnetic ions is close to a metal-insulator 
transition. From the theoretical point of view, the A- 
(BETS) 2 X and k-(BEDT-TTF) 2 X compounds have very 
similar band structures: in these three-quarter filled sys- 
tems, two bands are isolated from the two others by quite 
a large gap g. This can be interpreted as the sepa- 
ration between the bonding and anti-bonding orbitals 
on a dimer of molecules g. Projecting out the bond- 
ing orbital on each dimer, the system is thus effectively 
half-filled and therefore undergoes a metal-Mott insula- 
tor transition if U/t is increased M. Chemical pressure 
can change this ratio driving the system from a metal 
(or superconductor) to an insulator B. Replacing non- 
magnetic Ga 3+ by magnetic Fe 3+ , the electronic parame- 
ters change even less E3j. Even though this could also, in 
principle, drive the system from a metal to an insulator, 
this could not explain why the metallic phase is restored 
under magnetic field in a first-order transition. 

We now show that the magnetic character of the ions is 
important to drive the system into the insulating phase. 
Projecting out the bonding-orbitals from Ziman's model 
leads to a simpler 2-band model, with Hamiltonian: 

ij,er i i 

+ g a li B H Y^ S i + 9^bH Y, of 

i i 

where q creates a hole on the dimer at site i. Si is a spin- 
s' operator for the local moments. a\ = | ^ a o c ia a a f3Ci t 
(where a denotes the three Pauli matrices) is the spin- 
1/2 operator for the hole on site i. U and J are, respec- 
tively, the projected Hubbard repulsion and the Kondo 
coupling, iy is the tight-binding hopping integrals be- 
tween dimers [pj . g a and g are the g-factors of the local 
moments and itinerant electrons, respectively. 

Let us take the two limits of small and large U of this 
model: (i) At small U and J small enough, the phase 
is metallic due to imperfect nesting 0]. The localized 



spins are subject to an RKKY interaction. Treating 
the local moment spins classically, the total energy is 
E m etai — z X* ''x{Q ' af) S 2 , where z — 2 is the number of 
magnetic bonds, (ii) At large U, the system is a Mott- 
insulator. The electrons are antifcrromagnetically or- 
dered because of the Anderson superexchange process. 
Subsequently, the Kondo coupling forces the S=5/2 mo- 
ments to be antiferromagnetically ordered with respect 
to the localized electronic spins. The magnetic energy is 
— t;JS per site and the total energy of the AF Mott insu- 
lator (AFMI) is Eafmi — \JS. The gain in magnetic en- 
ergy is much larger in the Mott phase than in the metallic 
phase (J 2 x(Qaf) ~ J 2 jE F <C J, where E F is the Fermi 
energy). Let us assume that the expressions of the mag- 
netic energies are still valid for intermediate U |U| . If for 
J = 0, Eafmi > E meta i (the Ga compound is a metal) it 
is possible that Eafmi — \JS < E Metal — zJ 2 \S 2 , pro- 
vided that J is large enough or the difference between 
Eafmi and EMetai is small enough. A similar argument 
applies to the energy of the superconducting phase be- 
cause the RKKY interaction near Qaf is not modified 
in the superconducting state (l4| . 

Destruction of the insulating phase by temperature. 
Above the Neel ordering temperature (Tn ~ Jo) f° r the 
local moments the metallic phase has entropy of order 
ln(2S). In contrast, the insulating phase with AF order 
has zero entropy. Hence, to zeroth-order in Jo, the metal- 
insulator transition is first order and occurs at a temper- 
ature of T M i ~ (E meta i(J = 0) - E A fmi(J))/M 5 )- 

Destruction of the insulating phase by a magnetic 
field. We calculate the classical energies of the metal- 
lic and AFMI states as a function of the magnetic 
field. Doing this, we can neglect the electronic sus- 
ceptibility because Jo "C ti,4t 2 /U. (i) Metallic phase. 
We restrict ourselves to spiral ordering such as Si — 
(S cos a cos(Q. Ri), S cos a sm(Q.Ri),S sin a). The en- 
ergy is, E(H,a) = E Me tai- z-hS 2 cos2a- g a ^ B H S s\na. 
Minimizing this with respect to a gives: E(H) = 
E M etai - zJ S 2 - (g a fi B H) 2 /8zJ for H < H N = 
4zSJo/g a fiB which is the the critical field to align the 
spins, and E(H) = E Me tai + zJ S 2 - g a ^BHS for 
H > ij/v ■ (ii) Insulating phase. The energy is E(H, a) — 
Eafmi — ^JS cos a — g a jiBSH sin a. The minimization 
gives E{H) - Eafmi - \JS^l + (2g al i B H/J) 2 . Pro- 
vided that E^fetai + zJ S 2 < Eafmi, as the field in- 
creases the energy of the metal crosses that of the insu- 
lator, leading to a first-order transition into the metallic 
phase. 

Field-induced superconductivity. The argument for the 
Jaccarino-Pcter mechanism [RB is as follows. If the sys- 
tem is sufficiently two-dimensional, when a magnetic field 
is applied parallel to the layers, the orbital motion of the 
electrons is quenched. The upper critical field is then 
determined by the Pauli paramagnetic limit fig ]. If we 
first neglect the fluctuations of the localized spins and 



consider the regime where the moments are aligned by 
the magnetic field, the Kondo term in the Hamiltonian 
is replaced with JJ2i Si- ai — —JS^^ of. The effective 
magnetic field experienced by the electrons is H — H®, 
where H® = JS/{gjjbB) is a compensating magnetic field 
if J > 0. At H — H®, the Hamiltonian is the same as for 
the compound without the magnetic ions ( J = 0) at zero 
field. As A-(BETS)2GaCl4 is a superconductor, this map- 
ping shows that A-(BETS)2FeCl4 has to be a supercon- 
ductor as long as \H — H® \ < Hp, the Pauli limiting field. 
The nature of the superconductivity in the two materials 
should therefore be the same. This is supported exper- 
imentally by similar thermodynamic quantities in both 
compounds (T c Ga = 5.5 K and T c Fc = 4.2 K; ff£ a = 12 
T and H^ max - H e ~ 10 T). Tilting of the magnetic 
field out of plane giving a perpendicular component of 4 
T destroys the superconductivity Q], This value is com- 
parable to the upper critical field for A-(BETS)2GaCl4 
for perpendicular fields p). Note also that even if the 
magnetic field is in the plane, the orbital limiting field 
must be larger than H® to get superconductivity. This 
explanation gives J = 1.6 meV for H® = 33 T 0]. 

Effects of the fluctuations of the localized spins. 
The above argument neglects the spin flip terms 
J^ i (S' i + (Tr +S^cr^) in the Hamiltonian, where the +, — 
superscripts denote spin raising and lowering operators, 
respectively. Without the fluctuations, the two spin- 
states of the electrons have the same energy for H = H®. 
This is no longer the case when the spins fluctuate: the 
spin down can flip while the spin S z = — S is raised to 
1 — S at same time. Flipping of the spin up is, however, 
blocked because it would require lowering the spin of the 
S z = —S state. These processes renormalize the com- 
pensating magnetic field. To gain some insight on the 
relative importance of this effect, we consider the sim- 
ple problem of just one local moment and one electron. 
The compensating magnetic field is then given by (when 
g ~ g a ) H e — f§5^#e (this reduces to H® for small fluc- 
tuations, i.e. large S). The real value of J is therefore 
slightly larger than that extracted above. The second ef- 
fect of the fluctuations is to increase the on-site repulsion 
between electrons. Two electrons on the same site cost 
not only the energy U but also block the fluctuations be- 
cause the spin down is no longer allowed to flip. This 
extra repulsion is given by JS/ (AS — 2), which is negligi- 
ble compared to U . In summary, due to the large value 
of S, spin fluctuations associated with the local moments 
do not significantly change the physics. 

In order to more clearly establish that the field-induced 
superconductivity is due to the compensation effect, it 
is desirable to have an independent measurement of the 
exchange field. We now show how to extract H e from 
the Shubnikov-de Haas oscillations. In layered organic 
metals a magnetic field perpendicular to the layers will 
produce oscillations in the resistivity that can be related 
to the Fermi surface parameters |lq ]. In A-(BETS)2FeCl4 



at high magnetic field, the magnetic ions impose an ex- 
change field that splits the conduction bands (for spins 
up and down). We calculate the two corresponding fre- 
quencies that should appear in the oscillations |L7|. In 
the absence of an exchange field, as the magnetic field 
is tilted at an angle 9 away from the normal to the lay- 
ers, the oscillatory part is of the form cos[2tt F/(H cos 9)) 
where F is the oscillation frequency. The amplitude of 
the oscillations is proportional to the spin splitting fac- 
tor R s = cos(ttSo/2cos9), where the argument is pro- 
portional to the ratio of the Zeeman splitting to the Lan- 
dau level splitting, Sq — g*m* /m e , with renormalized 
mass and g- factor JL8| . In the presence of the exchange 
field, the spin-splitting factor is modified |l9[ . We get 
R s = cos[irS (H e /H- l)/2cos0]. The effect of this is to 
produce two oscillation frequencies, Fj cos9 ± SF where 
8F = S H e / '(4 cos 6). In A-(BETS) 2 FeCl 4 , Uji et al. ob- 
served two frequencies with a difference of 130 T/cos(9 
H . If we interpret the frequency difference as due to the 
exchange field [g0|, we extract H e — 32 T using the ob- 
served effective mass m* /m e = 4.1, and assuming g* = g 
pl[ . The interpretation of the splitting in terms of a 
corrugated three-dimensional Fermi-surface |9l is incon- 
sistent with the observed l/cos^ dependence because in 
that case the difference should vanish at the Yamaji an- 
gles 



LG 



22 1 , seen in the angular- magneto-resistance oscil- 
lations 1 9]7 Thus the magnetic oscillations imply that the 
compensating field should be about 32 T, in remarkable 
agreement with the optimal field for superconductivity. 

Electron spin resonance. It should be pointed out the 
frequency splitting discussed above occurs independently 
of the sign of J. However, the sign of J can be determined 
unambiguously by ESR. In the presence of the exchange 
field, the ESR frequency in the high-field regime, lo = 
gHB\H — H e \ P3l , will give H e and its sign. Note that, 
the magnitude of this frequency should be sufficiently 
small near H e that it is experimentally accessible. 

Based on the above picture and the analysis be- 
low we predict field-induced superconductivity in k- 
(BETS) 2 FeBr 4 . It is an AF metal below 2.5 K, and 
undergoes a superconducting transition at 1 K P4| , p5| . 
The magnetic oscillation spectrum also has two frequen- 
cies with a difference of 100T/ cos 9 and an effective mass 
of m*/m e = 8 Pq] . This gives an exchange field of 
\H e \ = 12 T. The critical field data for K-(BETS) 2 GaBr 4 
are not available; but we can estimate Hp from the crit- 
ical temperature assuming a BCS relation |L5| Hp a ~ 
1.8 kpT c / up = 1.2 T. With the above values for H e and 
Hp we would expect field-induced superconductivity in 
the range 11 to 13 T if J > 0. 

We now show how the upper critical field parallel 
to the layers can be greatly reduced when there is co- 
existing superconductivity and AF ordering of the mag- 
netic ions. This has been dramatically demonstrated in 
A-(BETS)2FeCLt under a pressure of 3.5 kbar. It is an 
AF metal above 3 kbar p7J] and undergoes a supercon- 



ducting transition at about 1 K pq ]. Normally, in lay- 
ered superconductors the upper critical field parallel to 
the layers is much larger than for the field perpendicu- 
lar to the layers. Here, the reverse happens! The upper 
critical field parallel to the layers is only H'j 2 = 0.05 T, 
whereas the perpendicular critical field is about 0.5 T 
p8| . This is in contrast with the Pauli limiting value 
estimated from the transition temperature, Hp = 2 T. 
We now show that this rapid destruction of supercon- 
ductivity by a magnetic field is due to the polarization 
of the magnetic ions and it can be related to the ex- 
change field. In the AF phase, the uniform component of 
the spins when a magnetic field is applied is (S Z )(H,T) 
leading to an exchange field: J(S Z )(H, T) . The staggered 
component of the spins superimposes an alternating mag- 
netic field that may open gaps in the Fermi-surface. As 
the nesting is imperfect, it will not destroy large parts 
of it and should not affect strongly the superconduct- 
ing properties. We will neglect this effect here. In this 
case, provided that the crystal structures of the com- 
pounds with and without the magnetic ions are similar, 
the upper critical fields of both compounds are related 
by \J{S z ){Hlf e {T),T) -g»BH!f e (T)\ = g» B H^{T). 
Measuring the upper critical fields and the magnetization 
curve allows a value for J to be extracted. For a classical 
antifcrromagnet with exchange Jo, the transverse magne- 
tization is given by g a fiBH/(4:zJo) at zero temperature. 

II Ff 

The relation then becomes |1 — g a /4zgJx(Q)\ H^ = 

II P 1 'p 

This shows that H^ can be much smaller than 
(because Jx(Q) ~ J/E F < 1). 

We now apply these ideas to K-(BETS)2FeBr 4 . The 
influence of the magnetic ions has previously been in- 
voked to explain why the upper critical field is anisotropic 
within the plane of the BETS molecules @. We 
rewrite the relation above between upper critical fields 
as |1 — H e /Hfq\ = Hp a /-ffj.2 e , having introduced the 
classical field to align the moments, Hn = AzSJo/g a (J-B 
p9J. This allows us to extract the parameter H e (or 
J) from the measurements of the critical fields. In k- 
(BETS) 2 FeBr 4 , H l }f e ~ f T for H || c, H N ~ 5 T || 
and Hp a ~ 1.2 T (see above). The positive solution is 
H e ~ 10 T, consistent with the estimate above. 

In conclusion, we have stressed the possibility of hav- 
ing a Mott insulator in A-(BETS)2FeCLi at zero magnetic 
field. The measurement of the charge gap as a function 
of field may help distinguish the Mott versus band insu- 
lator: for the Mott picture the gap should not vary signif- 
icantly with field whereas for the band picture it should. 
Furthermore, we have shown that the Hamiltonian that 
describes A-(BETS)2FeCU at high fields is simply related 
to that for A-(BETS)2GaCl4 with a compensating mag- 
netic field acting on the spins. We have interpreted the 
splitting of the magnetic oscillations as a signature of 
the exchange field, thus allowing us to extract the Kondo 
coupling. The strength of the exchange field equals that 
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of the optimal field at which superconductivity is ob- 
served. This strongly supports the Jaccarino-Peter effect 
and suggests that the nature of the superconductivity is 
the same in both materials. Note that this picture is 
true only if J > and can be confirmed by ESR. Using 
the same procedure, we have predicted that the organic 
superconductor K-(BETS)2FeBr4 should also exhibit a 
field-induced superconducting phase at about 10T, hav- 
ing extracted the Kondo coupling both from the magnetic 
oscillation splitting and from upper critical field data. 
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